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Recently published crater statistics on the small asteroids 25143 Itokawa and 433 Eros show a significant depletion of craters below 
approx. 100 m in diameter. Possible mechanisms that were brought up to explain this lack of craters were seismic crater erasure 
and self armoring of a coarse, boulder covered asteroid surface. While seismic shaking has been studied in this context, the concept 
of armoring lacks a deeper inspection and an experimental ground truth. We therefore present cratering experiments of glass bead 
projectiles impacting into granular glass bead targets, where the grain sizes of projectile and target are in a similar range. The 
impact velocities are in the range of 200 to 300 m s"' . We find that craters become fainter and irregular shaped as soon as the target 
grains are larger than the projectile sizes and that granular craters rarely form when the size ratio between projectile and target 
grain is around 1:10 or smaller. In that case, we observe a formation of a strength determined crater in the first struck target grain 
instead. We present a simple model based on the transfer of momentum from the projectile to this first target grain, which is capable 
to explain our results with only a single free parameter, which is moreover well determined by previous experiments. Based on 
estimates of typical projectile size and boulder size on Itokawa and Eros, given that our results are representative also for kms"' 
impact velocities, armoring should play an important role for their evolution. 

Keywords: Impact processes, Cratering, Regoliths, Asteroids, surfaces. Asteroid Itokawa, Asteroid Eros 
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1. Introduction 

Impact craters are ubiquitous on many solar system bod- 
ies and shaped their faces from their early evolution. Recently, 
images were taken of asteroids 25143 Itokawa visited by the 
Hayabusa spacecraft (Fujiwara et al., 2006) and 433 Eros vis- 
ited by NEAR Shoemaker (Veverka et al., 2000), providing 
an unprecedented quality as the spacecrafts closely approached 
and even landed on these small bodies. The images showed 
that while many other solar system bodies are saturated with 
craters, the surfaces of these two bodies are depleted in small 
craters. For Itokawa, the largest observed craters were 150 m 
across and the expected empirical saturation with craters (i.e., 
a surface coverage of around 25 %, Hartmann, 1984) was only 
found for the largest craters (Hirata et al., 2009). The surface 
coverage is declining for smaller crater sizes. On Eros, Chap- 
man et al. (2002) found an empirical saturation of craters larger 
than approx. 100 m but the same depletion of smaller craters 
as on Itokawa. Understanding the reason for the depletion of 
these small craters can be a clue to understand the physical pro- 
cesses on these asteroids and potentially even learn about their 
unknown interior. 

One hypothesis is that craters are being smoothed out by 
seismic activity, which might be driven by impacts. This has 
been extensively modeled by Richardson et al. (2005) with fo- 
cus on Eros. Laboratory experiments on this topic were started 
by Izenberg and Barnouin-Jha (2006) while more quantitative 
results are still needed. Based on these developments, O'Brien 
et al. (2006) were able to explain the depletion of small craters 



on Eros by seismic shaking and Michel et al. (2009) succeeded 
to explained it for Itokawa. However, some details remain un- 
clear and it is at this point not proven whether another eff'ect 
might not also play a vital role. Hirata et al. (2009) proposed 
that armoring of large boulders might explain the lack of craters 
on Itokawa. This effect, first introduced by Chapman et al. 
(2002, giving credit to F. Horz), builds on a surface covered 
with boulders, which are larger than the impactor size that would 
have been needed to form those craters. An impactor striking 
one of these boulders would then lose a considerable amount of 
its impact energy into the crushing or cratering of this boulder 
instead of forming a granular crater on the asteroid. A possi- 
ble direct evidence for this might be the bright spots on some 
of Itokawa's boulders which are believed to be small impact 
craters (Nakamura et al., 2008; Takeuchi et al., 2010). More- 
over, some boulders on the surface show evidence for being 
fragmented in place, i.e., on the asteroid surface (Nakamura 
et al., 2008). The physical processes and the efficiency of ar- 
moring are however widely unknown and subject of this exper- 
imental paper. 

Some experiments have been performed with relevance to 
our work. Most previous experiments of cratering in granular 
material involved projectiles which were larger or much larger 
than the target grain sizes. Cratering experiments with glass 
bead targets of varying grain sizes have been performed by Ya- 
mamoto et al. (2006). The projectiles in these experiments were 
much larger than the sub-millimeter target grains but a variation 
of the target grain sizes showed a variation of the crater size. 
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For this configuration, the effect can probably be explained by 
the cohesion of the grains, which had a diameter of 36 fim in 
the smallest case but nonetheless the experiment is noteworthy 
in our context because the coarsest grained targets (220 jum) 
will be directly comparable to our finest targets. 

An experiment where the target grain size is the same as 
the projectile is presented by Barnouin-Jha et al. (2005). They 
launched 3.5 mm glass bead projectiles at 1 - 2 km s"' into tar- 
gets of identical spheres. They found that the craters generated 
by 1 kms"' impacts were found to be very irregular in shape 
and only for the higher velocities had a conical shape. They 
also describe that small differences in the location of first con- 
tact between a projectile and target tend to determine the final 
shape of a crater An efficiency for possible armoring in these 
experiments is however not given. 

An example of an experiment to directly study armoring is 
given by Durda et al. (2011): granodiorite blocks of roughly 
5.75 cm edge length were placed on silica sand simulating as- 
teroidal regolith. The blocks were embedded to different depths 
from to 7.2 cm and exposed to a 3.2 mm aluminum projectile 
at 5 km s ' . For deeply embedded blocks (more than half of the 
block's edge length), the craters in the sand around the blocks 
were found to be of similar size as craters into pure sand (ac- 
cording to their Fig. 5) so they did not find significant armoring. 
Only for the slightly or not at all embedded blocks, the craters 
were about half of the size and featured many smaller craters 
produced by fragments of the granodiorite block. Durda et al. 
particularly studied the fragmentation of these blocks and found 
that the largest fragment after the collision was getting larger if 
the blocks were embedded to a deeper level, i.e., fragmentation 
of the block was less in that case. 

In spite of these experiments, there is no dedicated study to 
quantify the efficiency of armoring under a large variation of 
parameters. The open questions to be addressed are the follow- 
ing: at which size ratios (projectile to target grain) can armor- 
ing be efficient? What is the threshold size ratio for the set in of 
armoring? How efficient can armoring be for extreme size ra- 
tios? How does the morphology and shape of the impact craters 
change when armoring comes into effect? What is the physical 
reason for armoring? And finally, what contribution could ar- 
moring make on the observed crater statistics on small asteroids 
like Eros and Itokawa? 

To answer these questions, we will present experiments of 
glass bead projectiles impacting into monodisperse glass bead 
targets. The sizes are chosen so that the ratio between projec- 
tile size and target-grain size are centered around unity, which 
is referred to as t// by Housen and Holsapple (2011). In Sect. 
2 we will describe our experimental setup and the results are 
presented in Sect. 3. A discussion of the results, including 
a physical model for armoring under the conditions described 
here, will be made in Sect. 4. Moreover, we will directly apply 
our results to the asteroids Eros and Itokawa and sketch how it 
can affect and possibly explain the observed crater statistics. A 
final conclusion is drawn in Sect. 5. 



Table 1: Reference data of the 33 conducted cratering experiments, 
exp. projectile target-grain velocity crater crater 
no. diameter diameter diameter depth 

[mm] [mm] [m/s] [mm] [mm] 

1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 

14") 
15"' 
16"' 
17"' 
18"' 
19"' 

20 "' 

21 "' 
22 
23 

24 "' 

25 "> 

26 "' 
27"' 
28"' 

29 "' 

30 "' 

31 "' 
32 
33 



"' highly irregular crater 

"' strength crater in one grain or single crushed grain 

2. Experimental Setup 

We studied impact crater formation in granular glass bead 
targets within a reasonably wide range of target-grain to projec- 
tile-size ratio centered around unity. The target glass beads we 
used had diameters of 0.2, 1, 3, 10, and 30 mm and the glass 
beads used as projectiles had diameters of 1, 3, and 8 mm (see 
Table 1 for a listing of all experiments). All beads were from 
soda-lime glass and had densities of 2500 kgm"^. The 0.2 mm 
grains are the same as used by Yamamoto et al. (2006, TC, 220 

For the acceleration of the projectiles we used either a gas 
gun or a powder gun at Kobe University. The gas gun acceler- 
ates 3.2 mm glass beads to velocities of up to 290 ms"', where 
the impacts took place under normal air pressure. With a poly- 
carbonate sabot glued to the front of the 3 mm sphere, we were 
able to accelerate 1 mm glass beads to slightly lower velocities 
of up to 240 ms"'. With the powder gun we accelerated 8 mm 
diameter projectiles attached to a 15 mm (diameter and length) 
sabot to velocities in the same range, in this configuration under 
a reduced air pressure of 5000 Pa. In all cases the sabots were 
safely stopped before making contact to the targets. 
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Figure 1: Tiie velocity distribution of the impact experiments listed in Table 
1 . The experiments were performed in two series and the dashed lines corre- 
spond to the first series comprehending shots 3-17 while the solid lines show 
the distribution for the second series of shots 18-33. 




Figure 2: Example for the application of the spray-paint method in shot 27 with 
3 mm target grains. The crater is very irregular but the shape can clearly be 
recognized from the change in color. 



The projectile velocities were pointing downward so that 
the granular targets could be stored in steel bowls of 26 and 
33 cm diameter. To avoid any influence of the container, we 
chose the diameter such that it was always at least three times 
the expected crater diameter (the only exception are shots 32 
and 33). Only for the experiments 24-27 we used a 9 cm plastic 
bowl as the craters were very small in these cases. 

The projectile velocities were measured by a high-speed 
camera in back-light illumination, operated at 100,000 frames 
per second. A distribution of impact velocities, which are also 
given in Table 1, is shown in Fig. 1. As the velocities could not 
directly be measured in some experiments, we have to use mean 
velocities which are representative for these shots. A problem 
in doing so can be seen in Fig. 1: the experiments were per- 
formed in two runs with three month in between (dashed and 
dotted lines) and we found that the scatter in the velocities of 
the gas gun was small for each run (within 5 ms '), while the 



overall scatter is much larger (30 ms"'). The velocities were 
significantly slower in the second series. We would make a 
mistake by using overall mean values while the error from us- 
ing the mean values for the individual series is rather negligible. 
For the mean values of the experiments with the gas gun, we 
will use 229 and 287 for the first series, and 199 and 266 for the 
second (slower velocities for 1 mm projectiles). 

Shots 1 and 2 in Table 1 were neither considered in this dis- 
tribution nor for the calculation of the mean value. Those were 
the first shots in the first series (the experiment number also rep- 
resents the sequence order) and we attribute the lower velocity 
to dirt in the barrel. For the second run, we performed several 
test shots intended to clean the barrel and successfully avoided 
this scatter. The velocities in the powder-gun experiments were 
measured for each individual shot so that the significantly larger 
scatter (blue triangles in Fig. 1) is not a problem here. 

In the experiments performed with the gas gun, the for- 
mation of the craters were observed with a second high-speed 
camera with 5,000 frames per second, mounted at an inclina- 
tion of 45° with respect to the target. The crater diameters 
and depths were then measured by different techniques. (/) For 
most craters, we used a laser sheet (cp. Yamamoto et al., 2009) 
through the center of the crater and the 45° camera, which gives 
a precise value for the diameter (i.e., better than 5 %) as well 
as for the depth by simple trigonometry. This was however not 
possible for the craters formed in the 10 mm and 30 mm target- 
grain beds as well as for the small craters formed in the 3 mm 
target grains of shots 24-27. The reason is that in these cases, 
the size of the crater is in the order of the target grains and the 
shape is not clearly determined by the laser-sheet method. The 
laser is furthermore not simply scattered on the surface of the 
macroscopic particles but penetrates the target to be diffused. 
( ii) We therefore painted the upper layer of beads with a thin 
layer of spray paint (3 mm and 10 mm targets grains; see Fig. 
2) or marked each sphere of the upper layer with a marker pen 
(30 mm spheres; Fig. 4 in Sect. 3). From the movement of 
the first layer it was then possible to give a good estimate of the 
crater diameter, which can however still be a matter of debate. 
We therefore measured a minimum and a maximum crater size, 
denoted by a range of error bars (see Table 1 and Sect. 3). For 
the crater depth in these cases, we had to estimate the number 
of layers (mostly one), which gives a reasonable value within 
plus-minus one layer. (///) For the experiments 14 and 15, we 
measured the crater size with a ruler and found it to be consis- 
tent with but less precise than the spray-color method, (iv) For 
some experiments into the coarsest grained targets (shots 16- 
21 and 30-31) the formed crater was not a crater in the granular 
sense but rather a material-strength determined crater in the first 
target grain hit by the projectile. In these cases, the crater size 
in Table 1 denotes the crater in this grain. 

3. Results 

3.1. Crater morphology, diameter, and depth 

A qualitative overview of eight representative experiments 
is compiled in Fig. 3. The four pictures in the upper row show 
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Figure 3: Overview on how the crater size decreases with increasing target-grain diameter The upper row shows experiments with 3 mm projectiles while the lower 
row shows craters formed from the impacts of 1 mm projectiles. In the case of a 1 mm projectile and a 10 mm grain-size target, the crater is only a material-strength 
determined crater in the first grain. The scale is given by the size of the containers (26 and 33 cm diameter) and in the bottom right image by the diameter of the 10 
mm glass bead. 



experiments with 3 mm projectile beads and first of all it is ob- 
vious that the craters in these experiments are larger than those 
from the 1 mm projectiles in the lower row. The experiments 
with 8 mm projectiles are disregarded here as they do not fill 
the considered parameter range. Apart from the difference from 
the two projectile sizes, one can also see that the crater size is 
decreasing and getting fainter for increasing target-grain sizes 
from left to right. The impact crater formed by a 3 mm projec- 
tile into the 10 mm grained target (i/' = 0.3, top right) is hardly 
distinguishable and was therefore not considered in Table 1 (as 
mentioned above we then used the spray-paint method). On the 
other hand, an impact of a 1 mm projectile into the same target 
(i^ = 0.1, bottom right) does not produce an impact crater in the 
granular sense but rather leads to the cratering of a single grain. 
In shots 30 and 31, the first 30 mm target grain that was struck 
by the 8 mm projectile was completely fragmented but without 
a significant effect on the adjacent spheres - 0.3, Fig. 4). In 
these cases we define the volume of removed material as a crater 
although it does not resemble a known crater shape. It would 
moreover not be able to distinguish a crater like this under the 
natural conditions on an asteroid (without knowing the original 
surface). However, we take the crater diameter as one target- 
grain diameter and the depth of the depletion was observed to 
be approximately half of that. The image in Fig. 4 is a good 
example to visualize the armoring effect in our experiments. 

The so measured crater diameters and depths are plotted in 
Fig. 5 as a function of the target-grain diameter The slight dif- 
ferences in the impact velocities are not considered in this plot. 
As akeady observed in Fig. 3, crater size and depth increase 




Figure 4: In the impact of an 8 mm projectile onto a 30 mm target grain (shot 
31) the first grain that was hit by the projectile is largely fragmented but the 
adjacent grains hardly moved. The black crosses mark the top of the spheres 
before the impact. 



with increasing projectile diameter and decrease with increas- 
ing target-grain diameter. One can also see that the crater size in 
both plots remains relatively constant as long as the target-grain 
size is smaller than or equal to the projectile size. This is the ex- 
pected behavior from previous studies under conditions where 
we can ignore armoring. For some craters, the crater sizes are 
not so well determined as the craters do not have a clearly de- 
fined rim (see crater size measurement in Sect. 2), which is the 
case for granular craters in the coarse-grained targets (30, 10, 
and partly 3 mm target grains). In those cases, we estimated a 
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Figure 5: The crater diameter and crater depth for different target-grain sizes. 
For all three different projectiles it is clearly visible that the crater becomes 
smaller for large target-grain sizes while relatively unaffected as long as target 
grains are small enough. The filled symbols mark the strength crater where 
only one target grain was affected by the impact. 



minimum and a maximum crater size and this is denoted by the 
error bars. Our best estimate is the mean value of those two. 
For the experiments with fine-grained targets (0.2 and 1 mm) 
we could use the laser sheet and the error in the crater size is 
smaller than the natural scatter 

Also for the experiments with coarse-grained targets where 
we do not observe a granular crater but a strength-dominated 
crater inside a single grain (see Fig. 3 right bottom for an ex- 
ample), the crater size is well determined and no error bars are 
given. We precisely measured the diameter of these craters on 
photographic (30 mm grains) or microscopic (10 mm grains) 
images, and for the crater depth we made a geometrical assump- 
tion: we first computed the volume and the height of a spherical 
cap with a radius which is the crater radius. The crater shapes 
were either flat, in this case, the height of the spherical cap de- 
termines the crater depth. In another case, the crater was con- 
cave and could well be approximated by a second spherical cap 
mirrored at the flat surface. In this case, the crater depth is twice 
the height of the spherical cap. The according volume was later 
also used for the mass balance of the 10 mm grains. In shots 18 
and 19, the 30 mm target sphere did not only show a crater but 
also some mass loss on the opposite side. This was not consid- 



ered for the crater depth measurement but for the mass balance 
in the next paragraph. 

3.2. Non-dimensional scaling 

Qualitatively, the crater diameters in Fig. 5 all show a simi- 
lar behavior for the three projectile sizes. To make these results 
also quantitatively comparable, it would be desirable to apply a 
proper scaling to collapse all results into one curve. A standard 
method for the scaling of crater sizes is the non-dimensional n 
scaling, where instead of the crater radius the scaled radius 



1/3 

^« = I — I Re 



(1) 



is used (see, e.g., review by Holsapple, 1993), where nip is the 
projectile mass, is the crater radius, and pc is the density of 
the target material. We assume this density as 1600 kgm"^, 
which was measured by Yamamoto et al. (2006) for the 220 fim 
grained targets and is consistent with an expected random close 
packing porosity of 0.36. For smaller grains, Yamamoto et al. 
found a higher porosity, which can be accounted for by the on- 
set of cohesion. For the strength craters in the grains, we use 
the density of the bulk material (2500 kg m"^) for calculating 
the TTfi value. A scaled crater volume ny is also often used in 
this context and we also tried to plot this. However, due to the 
large scatter of some irregular shaped craters it did not lead to 
meaningful results here. The second important scaling parame- 
ter is the gravity-scaled size 



„2 



(2) 



which involves the gravitational acceleration g, the projectile 
radius Rp, and the impact velocity v. A standard way would 
now be to plot ttr against jt2, which then results in a power 
law with a characteristic slope and pre-factor for a given target 
material, e.g., glass beads. As these parameters do not include 
the target-grain size, we presume this power law ttr oc tt^" '^ 
(with the power 0.18 for glass-bead targets, Yamamoto et al.) 
and use the parameter ttk/tt^^^'^^, which should accordingly be a 
constant. 

These values are plotted in Fig. 6 and first of all we see that 
the three curves for three different projectiles collapse into one. 
We also see that the values are constant and roughly agree with 
the results of Yamamoto et al. (the final crater sizes as given 
in their Tables 2-4) as long as the projectiles are larger than the 
target grains {i// > 1; i/' is the projectile size normalized by the 
target-grain size). A small difference with the experiments of 
Yamamoto et al. is expected here as those were performed un- 
der reduced air pressure while our experiments are carried out 
under normal pressure conditions. Yada (2007) studied the ef- 
fect of air pressure on the size of the transient and final crater 
(same setup as Yamamoto et al., their TB target) and found that 
the final crater under reduced pressure below approx. 400 Pa 
is by 15-20 % larger than under normal pressure, which fully 
agrees with the difference between our data and the largest grain 
size for of Yamamoto et al. The difference should be signifi- 
cantly less for larger target grains and not play a role any more 
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Figure 6: According to the well-established tz scaling, the value n/dn^-^^ 
should be a constant. In our case, the value is constant for projectile diame- 
ters which are greater than the target-grain diameters but declines for smaller 
values. The black solid line represents the momentum transfer model described 
in the text. 



for the 10 or 30 mm grains. For target grains of similar size as 
the projectiles {if/ ^ 1), the mean value of the normalized crater 
size is already undoubtedly smaller and for even lower values 
(i^ < 1), the values steeply fall. 

For ij/ X 0.3, the values show a significant scatter, which is 
partly due to the fact that these crater sizes were the most un- 
certain compared to the other size measurement methods. But 
there was also a natural scatter observed, which for example 
manifests in the crater sizes of shots 25 and 26, where the dif- 
ference of the two crater sizes is as high as a factor of two. This 
was not a measurement error but already a stunning difference 
during the experiments, which is the reason that we reproduced 
it several times. This difference might be explained with the 
first point hit by the projectile: if the projectile hits more or less 
centrally on a target grain, the energy is released at a higher 
layer as when the projectile hits between the grains, penetrates 
into the granular medium and hits a target grain at a deeper 
layer. The deflected projectile fragments would have a greater 
influence on the crater formation. Part of this has also been 
described by Bamouin-Jha et al. (2005). 

It can moreover be observed that the scatter as well as the 
mean values of these t// » 0.3 data are increasing with decreas- 
ing absolute grain sizes. This can qualitatively be explained by 
a discrete effect due to the different target grains and their ra- 
tio to the crater sizes. In the case of the 30 mm target grains, 
the affected region was of the size of one target grain. If the 
target grain was centrally hit (which was the case), the next 
possible larger crater size would be three times bigger, because 
the adjacent grains would have to be displaced. In that case, 
a significantly higher energy would be needed and the n ratio 
would already be at the upper end of the scatter. This factor 
of grain size to crater size is decreasing with decreasing target 
grain size, as the crater size for the 10 mm target grains was 3 to 
4 spheres in diameter and the crater for the 3 mm target grains 



was 4 to 10 spheres in diameter. That means that a small scatter 
for the small target grains is easily possible while the scatter for 
the large target grains is less likely. While the point of impact 
as mentioned above is a reasonable explanation for the scatter 
in the shots with 3 mm projectiles, this effect should be taken 
into consideration for the larger target grains. 

It is important to note that the eight filled symbols encir- 
cled by the dashed line in Fig. 6 are those craters which are not 
due to a displacement of the granular medium but by cratering 
or fragmentation of a single target grain. In a strict sense, these 
values must not be plotted in this figure, which describes craters 
in the gravity regime. However, they can be interpreted as an 
approximated size of craters that might have formed on a gran- 
ular scale if this was possible. The upper size limit for those 
craters would be one target grain (when the adjacent grains did 
not move), the lower limit would be no crater. So we believe 
that this strength crater is a reasonable compromise for the size 
of a supposed granular crater. 

Although these craters are significantly smaller than the cra- 
ters in the fine targets, these values inhibit a surprisingly large 
scatter It is for example unexpected that the 1 mm projectiles 
(red circles) produce smaller craters in 10 mm glass beads than 
in 30 mm glass beads. As mentioned earUer, the impact param- 
eter was very small (central) in the case of the 30 mm target 
beads while it is unknown for the impacts into the 10 mm beads 
(i// - 0.1), which might explain this finding. The difference in 
the radius of curvature between the two grains would make a 
difference in the crater size that is opposite to what we find here 
(Fujiwara et al., 1993). 

3.3. Fragmentation of the first target grain 

We observed fragmentation of the target grains in most ex- 
periments with target grains larger than 3 mm. The smaller 
grains probably also fragmented but we did not see it as those 
were too many and we did not sift the targets after the impacts to 
find these. In the cases where we have a large granular crater, 
we were not able to find all fragments of these grains but in 
those shots where only one grain was affected by the impact, 
we either measured the crater size in this grain or the largest 
fragment resulting from the destruction from the grain. In Fig. 
7 we show the largest fragment (the cratered sphere is also a 
large fragment in a strict sense) as a function of the specific 
impact energy, i.e., the projectile energy normalized by the to- 
tal mass of projectile and the destructed target grain. As most 
of our collisions only produced a small crater, many values are 
close to unity. These data are shown in the inset plot at an 
enlarged scale. The black symbols show free collisions in the 
literature (Gault and Wedekind, 1969; Okamoto and Arakawa, 
2008), which will be further discussed in Sect. 4.5. For four of 
our shots with 30 mm target grains (shots 18-21), we were able 
to determine the normalized impact parameter (i.e., the offset 
from the central impact, normalized by the target-grain radius) 
by the slower high-speed camera mounted at 45°. These impact 
parameters were 0.15, 0.17, 0.07, and 0.17 and thus all reason- 
ably central. For the strength craters in the 10 mm grains (shots 
16 and 17), we were not able to measure the impact parameter. 
For the experiments where the 30 mm grains were disrupted 
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Figure 7: The mass of the first target grain's largest fragment (relative to the 
grain's original mass) as a function of the specific impact energy for eight indi- 
vidual shots (16-21 and 31-32). The values close to unity are also shown in the 
inset plot. A comparison to literature values (black symbols) will be drawn in 
Sect. 4. 



by the 8 mm projectile (shots 30 and 31), the impact parame- 
ter could not be determined due to a lack of the second camera 
but we tried to achieve a central collision, which is not a great 
challenge for 30 mm target grains. We thus also expect these 
impact parameters to be less than 0.2 as in the other cases and 
the data in Fig. 7 to be relevant for near-central collisions. 

4. Discussion 

4.1. Model description 

It is now desirable to understand the physical reason for the 
observed armoring and we will show that a surprisingly simple 
model based on momentum transfer is capable of explaining 
our results. If we consider the central collision between the 
projectile with mass and impact velocity v and the first target 
grain with mass nit, the target grain is not excited to the same 
velocity as the impact velocity if it is not much smaller than the 
projectile. From the conservation of momentum and energy, the 
velocity of the target grain would be 



(1 +e)v 
1 + mt/wp 



(3) 



if it was a free collision. Like in Newton's cradle experiment, 
this momentum would then be transferred to the subjacent grains 
and can thus be regarded as a particle velocity in the granular 
medium, which determines the pressure p = pcvg of the pres- 
sure wave (p and c are the density and sound speed of the gran- 
ular medium, see e.g., Melosh (1989)). The squared coefficient 
of restitution is a measure for the dissipation of energy and is 
unity for an elastic and zero for a completely plastic collision. 
Taking the relation between the n values, we can rewrite 



—0 1 8 



(4) 



where we assume that the velocity responsible for the crater for- 
mation is the particle velocity in the target and not the projectile 
velocity. This is not important as long as the projectile mass is 
much larger than the individual target grain mass (ni^ » Wt) as 
both velocities are in a constant relation in that case. However, 
if the projectile grain has the same size or is even considerably 
smaller than the individual target grains, this would become im- 
portant. In an independent approach to describe impacts into 
mortar-covered basalt blocks, Dohi et al. (2012) also found the 
particle velocity in the mortar as the controlling parameter to 
produce craters in the subjacent basalt rock (strength regime). 
Using Eq. (3), we can rewrite Eq. (4) as 



-0.18 



v2 



-0.18 



(1+fi) 

1 + ni^lnip 



0.36 



(5) 



where the first term on the RHS is the n2 value with the pre- 
viously expected projectile velocity and the second term is a 
new correction factor that describes the coupling between the 
projectile with mass ni^ and the granular target with grains of 
individual masses nii. The factor {\ + ef -^^ with e between 
and 1 varies only in a small range from 1 to 1.3. Therefore 
omitting this factor, we can write 



^-0.18 



1 + 



nit 



-0.36 



(6) 



This correction factor for the crater size is unity if the pro- 
jectile is much larger than the target grains and becomes ef- 
ficient in the armoring regime where this is not the case any 
more. With just one free parameter, the maximum value, Eq. 
(6) is plotted as a solid line in Fig. 6. The agreement is very 
good and only for values around tfr = 0.3 the theoretical curve 
systematically underestimates the crater size. This is however 
also a regime where we expect that specific conditions like the 
release point of the energy may be important (see above). The 
model is thus quantitatively capable to predict the relevant i/r 
value for the set-in as well as the width of the armoring regime. 

An approximation was made above that the projectile hits 
only one target grain, which is a good assumption if the projec- 
tile is small compared to the target grains. In contrast to that, 
if the projectile is larger than the target grains, it would interact 
with many projectile grains before being significantly slowed 
down. However, in that case, the velocity transferred to the 
target grains would anyway be the projectile velocity so that a 
further interaction would not make a difference for the initial 
velocity of the grain velocity Vg. 

4.2. Application to higher velocities 

An open question is still whether our results at velocities 
around 250 ms"' also hold for collisions at several kilometers 
per second as expected for the cratering on asteroids. These 
velocities are in the same range and often even higher than the 
sound velocity of glass or rock, while our velocities are signif- 
icantly smaller. One argument that our results are indeed also 



applicable to supersonic velocity impacts is the following: the 
velocities Vg of the target grains will be considerably smaller 
than the impact velocity v, for example by a factor of 1000 
which only needs a reasonable size ratio between projectile and 
target grain of 1: 10. Thus, even at a supersonic velocity impact, 
the target grains would only be accelerated to a fraction of the 
sound speed. 

On the other hand side, in the case of antipodal spallation of 
the first grain, the coupling to the second layer of target grains 
might be considerably different. This second layer would then 
not be accelerated by a target grain with its full mass but only 
by some smaller amount of spalled mass with higher velocity 
though. In any case, care should be taken as the momentum 
transfer at the grain boundaries might be very different for the 
hypervelocity collisions in general but this is so far unknown. 

We should also point out that a proper velocity scaling is 
not possible with our current experiments. It might for instance 
be conceivable that the horizontal axis in Fig. 6 includes some 
power of the velocity, which we would not see with the current 
data. In fact, a hypothetic scaling by the expected crater size 
Rc over the target grain size Rt gives a similar good fit for the 
three projectile size datasets as the one in Fig. 6. This would 



be proportional to R" 



i0.82,,0.36 



/Rt - the power of the target- grain 



size is close to unity and the velocity was nearly constant in our 
experiments. 

In this context, also a comparison with the experiments of 
Durda et al. (201 1) is interesting. These were performed at im- 
pact velocities of 5 kms"', which is only slightly below the 
sound speed of granodiorite (6 km s ' , Dayre and Giraud, 1986). 
They found armoring with a reduction of the crater diameter by 
a factor of two at most. Given the huge size ratio between the 
6 cm target blocks and the 3 mm projectiles, this is way off the 
range from our expectation from Eq. (6) - the craters are still 
too large. A velocity scaling in the horizontal axis of Fig. 6 
could solve that problem in principle but a power of 0.36 as de- 
scribed above would not make the data fit to our expectations. 
A plausible explanation of this misfit is the quality of coupling 
between the first target grain and the rest of the target. In the 
case of Durda et al., the granodiorite blocks were resting on the 
fine sand while our glass beads were supported by presumably 
only three contact points (see Sect. 4.5 for a deeper discussion 
on this coupling). Moreover, Durda et al. observed virtually no 
armoring when the blocks were deeply embedded into the sand 
and thus better coupled. 

A proper velocity scaling of the armoring effect is thus left 
for future work and the strong point we want to make in this pa- 
per is that we have to choose the particle velocity in the granular 
target Vg for the crater scaling instead of the projectile velocity 
V. The good agreement between our data and our momentum 
transfer model is still not a proof for the physical correctness of 
the model but nonetheless a good indication for this. 

4.3. Implication to crater statistics of Eros and Itokawa 

Assuming that the process works for high velocities as well, 
we can study the consequences from Eq. (6) in explaining the 
lack of craters on Eros and Itokawa as described in Sect. 1. 
First, we need to consider the boulder sizes on the surface of 



those bodies. For armoring to be efficient, these should be in 
the order of or larger than the impactors which would otherwise 
form the missing craters below 100 m. 

For Eros, Chapman et al. (2002) describe that boulders "< 
10 m in size dominate the landscape at high resolutions". More 
quantitatively, their distribution of boulders is presented in Fig. 
8 (gray open circles) in terms of the /?-plot. This value R is the 
surface coverage of boulders (usually used for craters. Crater 
Analysis Techniques Working Group, 1979) as presented in Fig. 
2 of Chapman et al. We consider here the boulders observed 
from a low altitude flyover while Chapman et al. also present 
a higher coverage for larger boulders in the images from higher 
orbits (i.e., lower resolution). Extrapolating these datasets, we 
arrive at a saturation of the surface {R - I) with boulders larger 
than around 2 m (Fig. 8) or 5 m (lower resolution data. Chap- 
man et al., 2002). This is indeed a similar size as an impactor 
which would otherwise form craters below 100 m as we will 
elaborate below. 

In the case of Itokawa we consider the boulder size distri- 
bution from Fig. 2 of Saito et al. (2006) and present it in the 
/?-plot in Fig. 8 (black open squares). For the error bars, we es- 
timated a 20 % error in the crater size from Fig. 2 in Saito et al. 
This distribution is significantly flatter than the distribution of 
boulders on Eros and if we extrapolate it, we arrive at a surface 
saturation with boulders of approximately 20 cm diameter with 
a large uncertainty due to the flat curve. Just as the different dis- 
tribution in the two datasets for Eros, we can also note here that 
the distribution of boulders on Itokawa depends on the exact lo- 
cation. There are very few boulders in the smooth MUSES-C 
or Sagamihara regions while many more can be found in the 
Ohsumi or Catalina regions. The local boulder sizes therefore 
differ from the global mean value given here but it should be 
sufficient for a rough estimate. 

Knowing a typical boulder size for the armoring of the sur- 
face, it is now interesting to see how Eq. (6) can help in describ- 
ing the crater statistics on these asteroids. We assume a flux 
of impactors following a differential distribution N(Rp) dRp oc 
Rp^ dRp that produce crater sizes according to 



1 + 



-2ff/3 



-a/3 



(7) 



written here in terms of projectile and target radius Rp and Ri, 
respectively, and K and a are parameters depending on the 
target material. In contrast to the parameters for glass-bead 
targets used before, we will here choose K - 0.69 and a = 
0.51 for sand or dry soil according to Schmidt and Housen 
(1987). We moreover assume an impact velocity in the main 
belt, where both asteroids are expected to have experienced 
their crater record, of 5 kms ' (Bottke et al., 1994a). If we 
use Eqs. (1) and (2) and consider the extreme cases with no ar- 
moring (Rt <K Rp) and distinctive armoring (Rt » Rp), we can 
simplify Eq. (7) to 

with p - \ - ajT) and fi - \ + 5q'/3, respectively. 
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Figure 8: Craters (filled symbols) and boulders (open symbols) on Eros and 
Itokawa (data taken from Chapman et al., 2002; Saito et al., 2006; Hirata et al., 
2009). The solid curve as partly described with Eq. (9) can represent the data 
when assuming a certain (constant) boulder size on the surface of the asteroids. 
The dashed line describes a boulder size which is a factor of two bigger or 
smaller than the best size. Caveats of this simple analytic model should be 
considered and are explained in the text. 



The surface fraction of craters (/?-plot) can be computed as 



(9) 



where we set the prefactor to match the empirical saturation that 
agrees to the craters of 100 m and larger on Eros. As a conse- 
quence of the new armoring term we find two different slopes 
(due to the different values of P) for two different regimes. The 
regime where we do not have armoring and (3 - l-a/3 shall be 
constant in the 7? -plot to represent the Eros data and we chose 
the only free parameter y = 2.66 to fulfill this constraint. While 
many authors assumed a slope in the range of y = 3.5 . . . 4.5 for 
this power based on dynamical modeling (Bottke et al., 1994b; 
O'Brien et al., 2006), new observations of asteroids of only few 
kilometers suggest a change in the power law and powers as 
small as 2.5-3 for these small bodies (Gladman et al., 2009; 
Ryan et al., 2012). In Fig. 8 we also show the crater statis- 
tics for the small craters on Eros (Chapman et al., 2002) and 
Itokawa (Hirata et al., 2009, all craters), in the same form as 
Fig. 1 (top) of Michel et al. (2009). The solid line is similar 
to Eq. (9) but representing the whole range without having to 
distinguish between the extreme cases. While this could not be 
explicitly solved, it is easy to tabulate with the computer 

With this easy approximation, we arrive at a slope of 1 . 1 in 
the armoring regime, which satisfactorily represents the data for 
Eros and the craters between 10 and 100 m of Itokawa. For both 
curves to match the data, we chose different (constant) boulder 
diameters of 5 m for Eros and 3 m for Itokawa. A change of a 
factor of two yields the two dashed lines in Fig. 8. This boulder 
size for Eros is comparable to the one estimated from the boul- 
der size distribution above while for Itokawa it is an order of 
magnitude larger. Seismic crater erasure might help to correct 
this but due to our very rough assumptions it is meaningless to 



force a better fit. We assumed here a surface covered with boul- 
ders all having the same size, while a distribution of boulders 
as presented in Fig. 8 would be more realistic and the boul- 
der coverage even depends on location. Both could naturally 
be studied with a Monte Carlo model while it is intricate with 
the simple analytic model presented here. Another caveat is the 
unknown applicability of Eq. (6) to the high-velocity regime 
(see Sect. 4.2 for a discussion). Therefore, a good fit in Fig. 8 
can only show up the potential of the armoring effect to explain 
the crater statistics on Eros and Itokawa while a more detailed 
study and further experiments in the high-velocity regime are 
still needed. 

4.4. Armoring on other solar system bodies 

Asteroid 243 Ida is another example for good quality im- 
ages as it was observed by the Galileo spacecraft in 1993. The 
smallest craters that could be resolved were of 100 m diam- 
eter and there was no depletion seen for any size (Chapman 
et al., 1996). Boulders were also found on Ida (Lee et al., 1996), 
which had sizes of 45-150 m while smaller boulders were not 
resolved. However, the surface coverage of the observed boul- 
ders was very low as there were in total only 17 boulders larger 
than 45 m. An extrapolation of the boulder distribution given 
by Lee et al. (1996) to estimate a size where it saturates the 
surface is not possible as the distribution is too flat. We there- 
fore cannot make a prediction whether Ida should be depleted 
in craters less than 100 m. 

To be briefly mentioned, there has also been a depletion of 
craters noted on our Earth Moon and on Mars' moon Phobos 
(Robinson et al., 2003). While this has been qualitatively ex- 
plained by seismic shaking, armoring would add another effect 
if the boulders are in the right size. On Phobos, 900 boulders 
have been identified which are greater than 5 m (Thomas et al., 
2000) and on our Moon it is generally known that the vicinities 
of large craters are populated with ejecta blocks, which may ac- 
count for a high density in some regions (Lee et al., 1986, and 
references therein). Without going into further details here, we 
want to note that armoring may give a contribution on many so- 
lar system bodies where we observe a depletion of small craters. 

4.5. Fragmentation of the first target grain 

In Sect. 3, we presented the amount of fragmentation of the 
first target grain in those cases where we had efficient armor- 
ing (Fig. 7) and concluded that the collisions can be regarded 
as quasi central. We will now draw a comparison to literature 
experiments on free collisions, which are denoted by the black 
symbols in Fig. 7. Gault and Wedekind (1969) presented su- 
personic velocity impacts into 7 to 10 cm spheres of decorative 
'crystal' glass. With a comparable density to soda-lime glass, 
the chemical composition given in the paper substantially dif- 
fers (e.g., high potassium content). Their data is shown by the 
black pluses. Okamoto and Arakawa (2008) performed impacts 
into soda-lime and quartz glass spheres of 8 to 13 mm at 1 to 
5 kms the results are shown by the black crosses. In their 
results, there is is no difference between the single soda-lime 
glass target and the other quartz targets so we do not distin- 
guish between those. The colored symbols are our collisions. 
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which were not free but it is expected that the grains are rest- 
ing on three subjacent grains, thus supported by three touching 
points. The Hterature values on the free collisions are in good 
agreement with our collisions of supported grains. 

The different datasets are not exactly comparable as the ma- 
terials were slightly different and also the velocities are not 
comparable. However, if we can conclude that it does not make 
a difference whether a grain fragments in free space or on the 
surface of an asteroid, this might have consequences for the in- 
terpretation of surface features on these bodies. We presented 
this comparison only to show our full set of data for others to 
use while only a dedicated study could show whether our sug- 
gestive conclusion is correct. 

The results of Durda et al. (2011) showed a significant dif- 
ference of the fragmentation strength of their granodiorite blocks 
depending on the level of embedding. This suggests that the 
quality of coupling to the secondary medium plays a role for 
the fragmentation and that three contact points are not an ad- 
equate coupling to show a significant difference to a free col- 
lision. To give a simple example with numbers, we imagine a 
2x2x2 cm-' glass block, resting on either three 10 mm glass 
beads or on 160,000 glass beads of 50 fim diameter The contact 
surface from the combination of the van-der-Waals attraction 
and the (earth) gravity of the glass block is given by Johnson 
et al. (1971, Eq. (19)). For typical material parameters of glass, 
we arrive at a total contact surface of 3.6 ■ 10"^ m^ for the 10 
mm beads and 3.2 ■ 10"** m^ for the fine 50 fim grains. The dif- 
ference is one order of magnitude and the reason is of course 
the huge number of grains, which are in contact with the glass 
block in case of the small grains. The coupling to the adjacent 
granular medium should thus also be significantly better in the 
case of Durda et al. so that this difference is not unexpected. 

5. Conclusion 

In the experiments at subsonic velocities presented in this 
article we showed that self armoring of the target plays a major 
role for the crater formation as soon as the target grains are in a 
comparable size range and bigger than the involved projectile. 
We presented a simple model based on the momentum transfer 
from the projectile to the first target grain, which is capable to 
describe our experimental results. This suggests that it is the 
particle velocity inside the granular medium excited by the pro- 
jectile and not the projectile velocity itself which is responsible 
for the crater formation. For projectiles which are significantly 
bigger than the constituent grains of the target, these velocities 
are in a constant relation and the effect cannot be observed. Our 
model has only one free parameter, which is the well known 
crater size for the situation where armoring does not play a role. 
Supporting DEM simulations would be interesting to test our 
model on a particle basis. These simulations could also easily 
study the effect of the impact spot, i.e., whether the projectile 
hits on top of one target grain or between the grains. This was 
found to make a large effect in the case where the projectile 
is in the size range of the target grains or slightly below (i.e., 
(// = 0.3..1, also see Barnouin-Jha et al., 2005). 



In describing the crater size distributions on small asteroids 
like Eros or Itokawa, it is still questionable whether our low- 
velocity results are applicable to high-velocity impacts on these 
bodies. Experiments in the high-velocity regime are needed 
here. However, if we can apply the results, it would be inter- 
esting to elaborate a full model (more advanced than the one in 
Sect. 4.3) including seismic shaking to reproduce the observed 
crater statistics and thus learn about the importance of seismic 
crater erasure compared to armoring. Another finding of our ex- 
periments is that craters in coarse granular media are irregular 
and hard to distinguish. This should also be taken into account 
when discussing the lack of small craters and might explain the 
even smaller fraction of craters smaller than 10 m on Itokawa, 
which is still overestimated in our approach. 

Finally, our results might have an impact on the Hayabusa 2 
mission, where it is planned to shoot a 10 cm projectile on the 
target asteroid. The current target of this mission is a C-type 
asteroid and although a previous observation of C-type aster- 
oid 253 Mathilde revealed no boulders on the surface (Thomas 
et al., 1999), the resolution of 160 m per pixel at its best might 
not have been sufficient to make a definite statement. In fact, 
possible positive mounds or block-like features have been de- 
scribed by Thomas et al. (1999) which were in the one or two 
pixel scale. The velocity for the impact in the Hayabusa 2 mis- 
sion will be 2 km s"', which is below sound speed, although we 
have to note that there are also some ordinary chondrites, for 
which the bulk sound velocity was measured below 2 kms"' 
(Yomogida and Matsui, 1983). 
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